Will all the data this has process we are able to plot information such as the change in net calcification between 2019 and the forcasted 2061
import pandas as pd
import numpy as np
import geoviews as gv
import hvplot.pandas
import geopandas as gpd
from bokeh.io import output_notebook
from bokeh.plotting import figure, show, output_file
# output_notebook()
gv.extension('bokeh')
df_all = pd.read_csv("Reef_Scatter.csv")
gdf = gpd.read_file("GBR Features Shapefile/Great_Barrier_Reef_Features.shp")
gdf
# only extract data that relates to reefs
reef = gdf['FEAT_NAME'] == "Reef"
reef = gdf[reef]
reef
sample = df_all.sample(n=9)[['Reef', 'Long', 'Lat', 'Sector','change (grams/m^2/day)']]
plot = gv.tile_sources.CartoLight * reef.hvplot(geo=True, color='gray',line_color=None) * gv.Points(sample, ['Long', 'Lat']).options(scaling_factor=20,
color_index='change (grams/m^2/day)', colorbar=True,
tools=['hover'], cmap='bwr_r', clim=(-3,3), size=10 , title="Change in calcification 2019-61"
).redim.range()
plot.opts(width=700, height=700)